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Abstract 

A well known behavior of EUV light curves of discrete coronal loops is that the peak 
intensities of cooler channels or spectral lines are reached at progressively later times. This time 
lag is understood to be the result of hot coronal loop plasma cooling through these lower 
respective temperatures. However, loops typically comprise only a minority of the total emission 
in active regions. Is this cooling pattern a common property of active region coronal plasma, or 
does it only occur in unique circumstances, locations, and times? The new SDO/AIA data 
provide a wonderful opportunity to answer this question systematically for an entire active 
region. We measure the time lag between pairs of SDO/AIA EUV channels using 24 hours of 
images of AR 1 1082 observed on 19 June 2010. We find that there is a time-lag signal consistent 
with cooling plasma, just as is usually found for loops, throughout the active region including the 
diffuse emission between loops for the entire 24 hour duration. The pattern persists consistently 
for all channel pairs and choice of window length within the 24 hour time period, giving us 
confidence that the plasma is cooling from temperatures of greater than 3 MK, and sometimes 
exceeding 7 MK, down to temperatures lower than ~ 0.8 MK. This suggests that the bulk of the 
emitting coronal plasma in this active region is not steady; rather, it is dynamic and constantly 
evolving. These measurements provide crucial constraints on any model which seeks to describe 
coronal heating. 
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1. Introduction 


A significant outstanding issue in current solar and astrophysical research is that of the 
heating of the solar corona. How is the corona heated to temperatures of greater than 1 MK when 
the photosphere below is only 6000 K? One observational approach to addressing this important 
question is to focus on particular areas in the corona such as active regions (AR). Often 
researchers have narrowed this problem even further by analyzing discrete coronal loops within 
an active region. A coronal loop is an intensity enhancement relative to the neighboring plasma, 
observed in extreme ultraviolet (EUV) or soft X-ray wavelengths. A standard analysis approach 
is to isolate a single loop in an image by subtracting contributions from background and 
foreground emission, and then examine the loop properties such as its differential emission 
measure distribution (DEM), temporal behavior, spatial structure, and densities (e.g. Schmeltz et 
al. 2001, 2011; Warren et al. 2002, 2008; Del Zanna & Mason 2003, Winebarger et al. 2003; 
Winebarger & Warren 2005; Ugarte-Urra et al. 2006, 2009; Tripathi et al. 2009; Reale 2010; 
Aschwanden & Boerner 2011; Viall & Klimchuk 2011). Such loop studies have revealed 
important information about their structure and heating. It is currently accepted by most 
researchers that coronal loops are composed of tens to hundreds of individual sub-resolution 
coronal strands, where a coronal strand is a miniflux tube for which the heating and plasma 
properties are uniform over the cross section (e.g. Cargill & Klimchuk 1997; Warren et al. 2002, 
2003; Winebarger et al. 2003; Winebarger and Warren 2005; Klimchuk 2006, 2009), although 
see Mok et al. (2008) for an different explanation. This important understanding has come about 
through the ability of hydrodynamic models of bundles of unresolved flux tubes to reproduce 
EUV and X-ray emission and light curves (e.g. Warren et al. 2002; Bradshaw and Cargill 2006, 
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2010; Klimchuk et al. 2008). In this framework, each coronal strand is heated by an energy 
impulse and then cools; many such nearby strands are heated within a finite time window (a 
nanoflare storm), creating a coronal loop. Importantly for the analysis we will present in this 
paper, the light curves of these loops exhibit a common time lag pattern: the intensity of a loop in 
a given EUV channel (or temperature band) will increase and then decrease, with the peak 
intensity of the light curves occurring in sequentially cooler channels as the coronal strands cool 
(e.g. Warren et al. 2002; Winebarger et al. 2003; Winebarger and Warren 2005; Ugarte-Urra et 
al. 2006, 2009; Viall & Klimchuk 2011). This is a common property of loops, even shorter loops 
observed in the cores of active regions (Viall & Klimchuk 2011). 

These analyses have contributed greatly to our understanding of coronal heating in loops, 
and are an important step towards understanding the heating of active regions and the corona in 
its entirety. However, it is always difficult to know how representative a given loop case study is. 
In fact, whether a group of heated strands is even identified as a coronal loop is likely biased by 
the temperature band of the observation and possibly the instrument resolution. Additionally, 
certain line of sight geometries preclude even the best background subtraction methods from 
fully isolating the emission from a single coronal loop (Aschwanden & Boemer 2011; Viall & 
Klimchuk 2011). 

Another limitation of these loop studies is that loops often make up only a minority of the 
total AR emission, with the diffuse emission in between the well defined loops contributing 60- 
90% of the total (Del Zanna & Mason 2003; Viall & Klimchuk 2011). Understanding this diffuse 
emission is therefore crucial for fully understanding the heating of active regions. It may be that 
the diffuse emission is simply composed of loops which are fading, or multiple overlapping 
loops, which are therefore indistinguishable. Indeed, Viall & Klimchuk (2011) presented a case 
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study of a location of diffuse emission in an active region whose light curves exhibit cooling 
patterns consistent with those of the previously mentioned loop light curves, suggesting a 
common heating mechanism. Alternatively, it may be that the diffuse emission is produced by 
coronal strands that evolve entirely independently, without the collective behavior that 
characterizes loops. The mechanism that heats these strands could be fundamentally different 
from the mechanism that heats the strands of a loop. For example, the strands which comprise 
the diffuse corona may be heated somewhat steadily and not undergo large heating and cooling 
cycles as loops do (Warren et al. 2010, 2011; Winebarger et al. 2011). Therefore, in this paper 
we expand on these previous case studies and analyze all of the coronal emission in an active 
region, rather than just single locations, in order to identify patterns of heat release and any 
subsequent cooling. 

Does the bulk of the coronal active region plasma exhibit cooling patterns as we have 
seen exhibited in countless loop studies, or is a different pattern observed for most of the active 
region? To address this important question, we measure intensity time lags between all possible 
pairs of six Atmospheric Imaging Assembly (AIA) EUV channels onboard Solar Dynamics 
Observatory (SDO) observed throughout an entire active region for 24 hours of data. Given the 
enonnous number of spatial pixels involved, it is not feasible to perform a comparison by eye as 
has been done for loop case studies. Additionally, it is desirable to develop a more statistically 
rigorous way to characterize the plasma evolution. Therefore we have developed an automated 
procedure that involves cross correlation of intensity time series. Upon applying this procedure 
we find ubiquitous signatures of plasma cooling in the bulk of the active region, consistent with 
the expected behavior of impulsively heated plasma. The measured time lags vary systematically 
as a function of location within the AR, which suggests that nano flares are more energetic in the 


4 



AR core. These measurements provide crucial constraints on any model of the coronal heating 
mechanisms operating in this active region, whether it is nanoflares, or something else. 


2. Methods and Analysis 

We analyzed 24 hours of EUV data for NOAA active region 1 1082 taken with SDO/AIA 
(Boemer et al. 2011; Lemen et al. 2011). We plot the normalized temperature responses 
functions (adapted from Viall & Klimchuk 2011) of the 6 SDO/AIA EUV channels in Figure 1 : 
we show 131 (black), 171 (cyan), 193 (orange), 211 (blue), 335 (green) and 94 (red). The 94 and 
131 channels are bimodal in their temperature sensitivity, with peaks at both low and high 
temperatures, while 335 is broad and 171, 193 and 211 are singly peaked and narrow. These are 
the published response functions at the time of this writing; however there is evidence that the 94 
channel in particular may have additional low temperature contributions (e.g. O'Dwyer et al. 
2010; Aschwanden & Boerner 2011; Foster & Testa 2011; Testa et al. 2011). These features of 
the response functions (e.g. bimodality) are important for interpreting our results. As shown in 
Viall & Klimchuk (2011), as an impulsively heated loop strand cools, we expect the intensity to 
peak in the different AIA channels with different orderings that depend on the nanoflare energy. 
Specifically, we found that for weaker nanoflares, we expect the cool contributions to the 94 
channel to dominate the light curves, and 94 should peak after 335, 21 1 and 193, but before 171 
and 131. For strong nano flares where the plasma is heated to much higher initial temperatures, 
we found that the hotter component of 94 will dominate, in which case we expect 94 to peak 
before all of the other channels. Regardless of the heating mechanism, if the 94 light curve rises 
and then falls preceding all of the others, the hot component dominates, however if the 94 light 
curve peaks between 193 and 171, then the cooler component dominates. 
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It is important to note that for impulsive heating, it is the emission from the cooling 
strands which dominates the light curves (Bradshaw & Klimchuk 2011). There is relatively little 
contribution from strands as they undergo heating because the heating phase is short lived and 
because the densities are relatively low. By the time chromospheric evaporation fills the strands 
with dense plasma making them visible, cooling is already well underway. 

The active region we analyze here was located in the northern hemisphere, near disc 
center on 19 June 2010. In Figure 2 we show images taken at 3.5 UT in the 6 EUV AIA channels 
described in Figure 1. From left to right, top to bottom we display the 94, 335, 211, 193, 171 and 
131 A channels. The field of view is 450x450 pixels, where 1 pixel = 0.6”. The image intensities 
are all on a linear scale. This active region has a bipolar structure, was not associated with a 
sunspot and did not have any flaring activity above B -class for the time period that we analyzed. 

We use the 30 second cadence, level 1.5 data for our analysis, for which all of the 
channels are coaligned, and have the same plate scale. We derotate the data so that the active 
region emission lies in the same pixels throughout the 24 hour period. To do this, we make the 
simplifying assumption that the whole active region rotates at one rate, and that the rotation can 
be approximated by the linear projection of the angular rate in the plane of the sky. These are 
reasonable since the AR is near disc center and has a small latitudinal extent. Any distortion of 
the coronal structures due to differential rotation originates at the magnetic footpoints, and these 
footpoints are all rooted within 10 degrees latitude of each other. While we do not claim to have 
completely eliminated the effects of differential rotation, a more involved correction is not 
worthwhile, since the active region will experience other forms of evolution during a 24 hour 
period. 
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For the latitude of the middle of this active region, the Snodgrass empirical rotation rate 
(Snodgrass, 1983) is equivalent to 300 seconds per AIA pixel. With this rotation rate, we align 
all images with the first image on 19 June 2010; additionally, we account for fractional pixel 
rotation, in which emission located entirely in one pixel in the first image is split between two 
neighboring pixels in subsequent images. For each image we use the time change since the first 
image to compute the fractional number of AIA pixels that the AR has rotated. For example, the 
12 th image is taken 330 seconds after the first image, and the AR has rotated 1.1 AIA pixels. We 
move the entire image over the integral number of pixels (for this example, 1 pixel), and we 
linearly interpolate between that pixel and its neighboring pixel to estimate the intensity change 
due to the remaining fractional pixel rotation. Linear interpolation to find the additional intensity 
change due to the fractional pixel rotation is crucial for time series analysis such as we perfonn 
in this paper. If the images are aligned only by integral pixel amounts, a periodic intensity 
change will be introduced into the time series equal to the rotation rate in seconds per pixel (for 
this data set, 300 seconds/ AIA pixels). Note that any slight errors in our estimation of the 
rotation rate or differential rotation effects, will take place over time scales much longer than the 
transient features important for this analysis, and will not affect the results presented here. 

At each pixel, in each channel, we construct a 24-hour time series from these coaligned, 
derotated images, which we then subdivide into smaller time windows (two consecutive 12-hour 
time series and twelve consecutive 2-hour time series). For each of the 12-hr intervals and 2-hour 
intervals, we construct time-lag maps in the following manner. We use the IDL procedure 
(/CORRELATE. PRO to compute the cross correlation value of the intensity time series at one 
pixel in one channel with the same pixel time series observed in another channel. We compute 
the cross correlation values as a function of temporal offsets (up to 1 hour for the 2-hour time 
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series, and up to 2 hours for the 12 hour time series) for both positive and negative offsets. We 
identify the time offset at which the peak cross correlation value is reached as the time lag, and 
construct time lag maps with the time lag measurements. 

In Figures 3 and 4 we illustrate an example of this technique applied to one pixel in this 
AR. In Figure 3 we show the normalized light curves from 2-4 UT (7200- 14400s) observed at 
the location indicated with the arrow in Figure 2. From bottom to top we plot the light curves 
observed in channels 131, 171, 193, 211, 335, and 94, each offset by 0.5 in y. There is at least 
one full heating and cooling cycle on strands along the line of sight visible in this time series. 

The intensity rises and falls in all of the channels, with each light curve temporally offset as a 
function of temperature. The peak intensity is reached in 335 around 9000s, 211 peaks around 
10000s, 193 peaks around 10300s UT, 94 peaks around 10600s, and 171 and 131 peak around 
1 1000s. This dynamic, time lag behavior is consistent with previous analyses of loop light curves 
as we described in the Introduction. Both the beginning and the end of this time series appear to 
capture cooling plasma as well. At the beginning of the time series, 211, 193 and 94 all decrease 
in intensity relative to their initial intensities, and 171 reaches peak intensity shortly after the 
beginning of the time series, and then also fades. This suggests that this time series began after 
the plasma has cooled to ~2 MK from an unknown maximum temperature at an earlier time. 
Similarly at the end of the time series, all of the channels exhibit intensity changes with time lags 
between hotter and cooler channels that are clearly associated with cooling plasma; however the 
time series ends before the intensities in channels 94, 171 and 131 have time to fade fully. 

In Figure 4, we show cross correlation values as a function of offset time between pairs 
of the light curves shown in Figure 3. The offset time is the amount of time the second light 
curve is shifted relative to the first; a positive time lag indicates that the second light curve 
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followed the first light curve. For cross correlating the channel pairs, we generally put the hotter 
channel first (though 94 and 131 could be both hotter and cooler than other channels due to the 
bimodal nature of their response functions). The choice of which channel is first is trivial, as the 
opposite choice will give the same time-lag value, but of opposite sign. The peak cross 
correlation value for 21 1-193 (green) is 0.9 and is reached when the 193 light curve follows the 
211 light curve by 300s (5 minutes). We record +300s at this pixel location for this channel pair 
in our 21 1-193 time lag map. The 335-211 (blue), 335-193 (orange), and 335-171 (cyan) cross 
correlations show similar patterns, with the peak cross correlation value reached when 21 1, 193 
and 171 each follow the 335 light curve at 700, 1000, and 1500s time lags, respectively. This is 
exactly what we would expect for a cooling plasma, given the respective temperature 
sensitivities of these channels and the approximate delays between the time series (Figure 3) that 
we identified by eye. 171-131 (black) reaches peak cross correlation value at no time offset 
between the two channels, and we record zero for the 171-131 time lag map. Lastly, we show 94- 
335 (red), which attains its peak cross correlation value at negative 1 100s, indicating that 335 
actually precedes 94, rather than follows it. This is consistent with the behavior we identified by 
eye in the time series. In the context of a nanoflare storm, this is consistent with a weaker heating 
scenario, as shown in Viall & Klimchuk (2011), in which the cooler component of the 94 
channel dominates. We include the entire 2 hour time series for this cross correlation example, 
just as we do when we apply this technique to the entire data set, even though the beginning and 
end of the time series may not coincide with the beginning and end of the plasma dynamics. 
Though this time series seems to begin while plasma is already in the middle of cooling, the light 
curves exhibit the same hotter-to-cooler time lag behavior throughout the time series, therefore 
the time lag signature comes through in our analysis of the cross correlation values. 
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3. Results 


We perform the analysis described above to every pixel in the image set, regardless of 
which type of AR structure is present (e.g., fan loop, AR core, long loops, short loops, hot loops, 
warm loops, or loop foot points). We compute the time lags between all possible pairs of the 6 
EUV channels for all of the 12-hour and 2-hour datasets. We display the results in the form of 
time lag maps where each pixel value is the time lag associated with that particular pair of 
channels over that time window. For all of the maps we show a color bar, which indicates the 
time lag range. Blues, greens and blacks indicate negative time offsets where the second channel 
precedes the first channel; reds, oranges and yellows indicate positive time offsets, where the 
second channel follows the first channel; the olive green color indicates that within the data 
resolution there is zero temporal offset between the two channels. 

3.1. The 12-hr time series 

We begin by presenting the results from analyzing the first 12-hour dataset. We display 
maps for every channel pair in Figures 5a and 5b. At the top of each panel we indicate the 
channel pair used . All panels are on the same color bar, except 211-193 and 171-131, which have 
a steeper gradient to highlight the small, but nonzero, time lag between those channels. In all 15 
panels, a clear time lag signal persists. Notice that though we unifonnly apply the same analysis 
to all of the pixels in the image, the structure of the AR is apparent in all of the maps. In the 335- 
211 panel (right column, middle row, Figure 5a) the active region is dominated by positive time 
lags, with some zero time lag, but very few negative time lag pixels. The pervasiveness of 
positive time lags indicates that throughout most of the active region the 335 light curve 
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variability precedes the 211 variability; only very rarely does the 211 channel precede the 335. 
This is consistent with a scenario in which the majority of the coronal active region plasma is 
heated to temperatures of at least 3 MK and cools to temperatures lower than 1 .6 MK. Recall that 
for impulsive energy release, the heating phase is much fainter than the cooling phase and does 
not influence the light curves. The other three 335 maps (bottom row Figure 5a: 335-193, 335- 
171, and 335-131) follow a pattern consistent with that of the 335-211 map, namely positive time 
lags dominate the active region, largely in the same regions where they dominated the 335-211 
map, with some zero time lag pixels, and very few negative time lag pixels. Comparing with the 
images of the AR in Figure 2 it is clear that the zero time lag locations are mostly coincident 
with moss, which is the transition region footpoints of hot coronal strands. The other location 
with many zero time lag pixels occurs in association with fan loops, though the pattern is not 
ubiquitous across all of the channel pairs or even throughout all of the fan-loop region. Taken 
together, these 4 maps indicate that throughout most of the active region, for these 12 hours, the 
335 light curves precede the 211, 193, 171 and 131 light curves; plasma that reaches 
temperatures hotter than 3 MK is cooling down to well below 1 MK. In general, the time lags are 
shorter in the 335-211 map, longer in the 335-193 map, and longest in the 335-171 map. This is 
expected for a cooling plasma, as 21 1 and 193 have hotter peak sensitivities than 171. 

In the top and middle row of Figure 5a we display the 94 pairs: 94-335, 94-211, 94-193, 
94-171 and 94-131. As we discussed earlier, the 94 channel has significant sensitivity both at 7 
MK and 1 MK. Nevertheless, we always compute time lags with 94 as the first time series in the 
pair. The 94-335 map mostly exhibits negative time lags in the active region, indicating that the 
335 light curve precedes the 94 light curve, just as in the example pixel of Figures 3 and 4. In 
contrast, the core of the active region has predominately positive time lags, indicating that the 
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335 light curves follow the 94 light curves there. This suggests strong nano flares in the core and 
weaker nano flares outside. We see a similar pattern for the 94-211 and 94-193 pairs of positive 
time lags observed outside of the core of the active region, with negative time lags observed in 
the core. The negative time lags for the 94-211 map and the 94-193 map are smaller than 94-335, 
as expected from their cooler peak sensitivities. The 94-171 and 94-131 maps differ from the 
other three: in these maps the vast majority of the active region is dominated by positive time 
lags, with few negative time lag pixels. In all of the 94 maps the time lag values at the loop 
footpoints, or the moss, generally are consistent with zero time lag. 

In Figure 5b we show the 6 remaining channel pairs: 211-193, 211-171, 211-131, 193- 
171, 193-131 and 171-131. The 211-193 map exhibits almost exclusively a positive or zero time- 
lag, with 193 following 211, and there is a clear tendency for the time lag to increase with 
distance from the core. Importantly, there is virtually no occurrence of negative time lags (193 
preceding 21 1), even very small ones. The 211-171, 211-131, 193-171 and 193-131 maps show 
similar patterns to those present in the 335 maps, namely the cooler channels follow the hotter 
channels (indicated with a preponderance of positive time lags) for the majority of the active 
region. In the moss regions we find time lags consistent with zero in these maps too. Lastly, in 
the 171-131 time-lag map (also with a compressed color table) almost the entire active region is 
correlated near zero time lag. There are some regions with slight positive time lag, and there are 
some regions with slight negative time lag, but zero time lag pixels dominate much more than in 
any of the other maps. 

We repeat this analysis on the second 12 hours on 19 June 2010 and find qualitatively 
similar results in the maps: positive time lags dominating ah but the first three 94 maps; 
progressively larger time lags as the peak temperature sensitivities of the channel pairs are 
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further apart; a 94 positive-to-negative inversion signature in the core of the active region; and 
zero time lag at the moss footpoints. 


3.2. The 2-hr time series 

We repeat the analysis, dividing the data into twelve consecutive 2-hour windows. This 
tests the degree to which single, short-lived intensity variations, surrounded by otherwise steady 
emission, dominate the 12 hour results. Do the 2-hour results display qualitatively the same 
results? We display time lag maps in Figure 6 computed using the first 2 hours of data in the 
same format as the 12-hour maps. We calculate cross correlation values for up to an hour 
temporal offset in the positive and negative direction, rather than a full 2 hours, for these shorter 
time windows. As the temporal offset increases, the number of data points used to compute the 
cross correlation decreases. Any temporal offsets greater than half the data window (in this case, 
greater than 1 hour), are computed using the variation of less than half of the time series, and 
therefore much less meaningful. 

Qualitatively these results are the same as the 12 hour results, though there is more noise 
in these results. The 94-335, 211 and 193 channel pairs exhibit negative time lags in most of the 
active region, with positive time lags in the core, while the other 6 panels in Figure 6a are 
dominated by positive time lags throughout the active region. In 6b we see that positive time lags 
persist for the majority of the active region in the 21 1-171, 21 1-131, 193-171 and 193-131 maps. 
Again, the 2 1 1 - 1 93 and 171-131 maps (Figure 6b) are displayed using a steeper color gradient to 
highlight the small but non-zero time lags. In 21 1-193 there is a larger occurrence of negative 
time lags than in the 12-hr maps, but the map is still dominated by positive time lags. The 171 - 
131 map has fewer olive green pixels than the 12-hour version, but as in the 12-hour version 
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there is a mix of positive and negative time lags. All of the maps have zero time-lag correlations 
in the moss, as in the 12-hr maps. 

These patterns persist for all 12 of the 2-hour maps in this 24 hour period. To illustrate 
this, in Figure 7 we show the 2-hour 335-211 maps for all twelve intervals in our dataset. All 12 
maps are consistent with a pattern of ubiquitous cooling for the duration of the 24 hour period. 
The details change across the different panels, and the exact time lag value at a given pixel may 
vary from map to map. However, a basic pattern of positive time lags persists across the entire 
24 hour period. In general, pixels for which positive time lags are measured have positive time 
lags in all 12 of the maps and regions of zero time-lag almost exclusively occur near the moss. 

There is a trade-off between window length and noise suppression. Comparing Figure 5 
with Figures 6 and 7 it is clear that the 2-hour results contain more noise, simply due to having 
fewer total points per time series to cross correlate. Additionally, the 2 hour windows are more 
susceptible to ‘edge effects’, where the time series captures only part of the plasma variability. 
This occurs because the start time of the windows are arbitrary (every two hours, on the hour), 
and may have nothing to do with the physical beginning of any dynamics. For example, in the 
nanoflare storm scenario, if a pixel time series contains only the very end of one nanoflare stonn, 
followed by only the very beginning of another, spurious time lags may result. In such a pixel the 
brightening and fading of 21 1 may even ‘precede’ that of 335, simply because the 335 peak from 
the first stonn occurs before the time window begins and is missed, and likewise the 211 peak 
associated with the second stonn occurs after the time window ends and is missed. As the time 
window increases, the likelihood of capturing many full cycles of plasma variability increases, so 
these spurious edge effects will be less and less significant. 
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The 12 hour maps have much better signal to noise and reduced edge effects, however 
there may be a question of persistence of the cooling. Perhaps the 12-hour maps are dominated 
by the time lag cooling signal simply because at some brief time in those 12 hours there was a 
feature which cooled through those pixels. With our 2-hour maps we can immediately rule out 
this scenario. If this were the case, then we would expect only the 2-hour windows that contained 
the transient feature to exhibit the time lag cooling signature. The 2-hour window results (Figure 
7) demonstrate that it is generally not a single, brief, transient event (such as one very bright 
loop) along the line of sight which dominates the full 12 hours. We see qualitatively similar 
behavior for all of the sub-windows for the entire 24 hour time period demonstrating that the 
cooling behavior is continual. 


4.Discussion 

Our main interpretation of these measured time-lags is that the plasma must be cooling 
throughout most of the active region for the full duration of the 24 hours. Furthermore, because 
positive time lags are widespread for the 335-21 1, 335-193 and 335-171 time lag maps with the 
longest time lags observed in the 335-171 maps, it is clear that at least some plasma must be 
heated to temperatures of greater than 3 MK, and must be cooling down to temperatures less than 
0.8 MK before undergoing any further heating. If there were positive time lags only present in 
the 335-211 and 193-171 maps, but not in the 335-171 maps, that could indicate two separate 
populations of plasma: one population which cools down from temperatures greater than 3 MK, 
but which heats back up before it cools below 1 .6 MK, and a different population which never 
gets hot enough to emit strongly above 1 .6 MK before cooling down. 
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These observations of the entire active region are consistent with the known cooling 
behavior of coronal loops, as we discussed in the Introduction, as well as consistent with our 
nanoflare storm models coupled with the AIA response functions discussed here and shown in 
Viall & Klimchuk (201 1). They are also consistent with a recent study of the distributions of soft 
X-ray intensity fluctuations (Terzo et al. 2011) which suggests widespread cooling. Even the 
more complex behavior exhibited by the 94 pairs is easily understood in the context of the 
nanoflare models, the bimodal nature of the 94 temperature response, and the other, more 
straightforward maps. As we discuss in Viall & Klimchuk (2011), the only way for 335, 211 and 
193 to precede 94 (the negative time lags in those maps) but for 171 and 131 to follow 94 
(positive time lags) is if the cooler, 1 MK component of 94 is dominating the light curves of 
those pixels. Likewise, the only way for 335, 211, 193 and 171 to all follow 94 is if the hot 
component of 94 dominates the light curves of those pixels. All of this suggests that the plasma 
in the longer loops, away from the core of the active region have 94 light curves dominated by 
the 1 MK plasma, while the shorter, core loop light curves tend to be dominated more by the 
hotter component of 94. Recalling the channel order predictions of Viall & Klimchuk (2011), 
this implies stronger nano flares in the core of the active region and weaker nanoflares outside. 
This is not surprising, since nanoflares are likely to be magnetically driven, and the magnetic 
field strength decreases with distance from the center of the active region. In contrast to the 94 
channel, we only seem to find evidence for the cool component of 13 1, suggesting that the 
coronal plasma in this active region does not strongly emit at temperatures greater than ~ 10 MK. 
Such ultra-hot plasma either has a very small emission measure or is strongly influenced by 
nonequilibrium ionization effects (e.g., Reale & Orlando 2008; Bradshaw & Klimchuk 2011). 
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The nanoflare model naturally explains the observations presented in this paper, 
including the time lag inversion in the 94 maps. What about other heating scenarios - can they 
too reproduce all of the aspects of the time-lag maps? For example, how could steady heating in 
which the plasma is constantly maintained at a single temperature and never has a chance to cool 
produce such time lag signatures? In this case the emission is by definition time stationary, and 
therefore the variability and resulting cross correlation will be exclusively due to noise. Noise 
will produce approximately equal amounts of positive and negative time lags, those time lag 
values will be random, and zero time lag will occur no more often than any other particular time 
lag value. Furthermore, our maps exhibit clear spatial organization in that the basic AR structure 
is visible; noise is not expected to have any particular spatial organization. Lastly, time lags 
resulting from noise would not produce the observed increase in time lags as a function of 
channel separation. We can be confident that for this active region, for these 24 hours, the 
majority of the emission cannot be due to truly steady heating. 

There are other heating scenarios other than impulsive nanoflare heating which may also 
be able to reproduce some of the time lag behavior. For example, truly steady heating can 
produce time-variable emission under thermal non-equilibrium conditions (e.g., Klimchuk, 
Karpen & Antiochos 2010; Lionello et al. 2011). Another possibility, type II spicules, which 
inject heated chromospheric material into the corona (De Pontieu et al. 2011), also contribute hot 
plasma that may produce a cooling signature such as those we observed here. A third possibility 
is a variation on the basic nano flare model, in which the heating is quasi-steady in the sense that 
the nanoflares repeat frequently enough on each strand so that the plasma only cools partially 
before being reheated. This would not, however, explain the observed time lags between 
channels that are widely separated in temperature. In our companion paper we model a range of 
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nanoflare frequencies for direct comparison with these time lag results. In that paper we also 
consider the possibility of nanoflares initiating randomly on physically separate coronal strands 
along a single line of sight, without the coherence expected for a bundle of strands comprising a 
loop. Regardless of the outcome, the results we present here are important constraints on all 
models of heating of the active region coronal plasma. 

Finally, there are significant moss/footpoints areas where the time lags are consistent 
with zero in most if not all of the channel pairs. This result is expected for the transition region 
emission from impulsively heated strands. For individual strands, the light curves in a particular 
channel from cooling coronal plasma are narrow, and the peaks have a clear separation in the 
different channels. This is not true for the emission coming from the transition region. Transition 
region light curves are much broader, and because the intensity tends to scale with pressure at all 
temperatures throughout the transition region, the channels all peak at roughly the same time 
(when the pressure peaks). This results in much greater overlap and temporal correspondence of 
all of the light curves. Consequently, we expect time lag signatures to be much less pronounced 
in the footpoints of impulsively heated strands than in pure coronal observations, just as we 
observe in our time lag maps. The other location in the AR that exhibited time lags consistent 
with zero was the fan loop area, though it was not as pronounced as the moss area. A possible 
cause is transverse waves carried on those strands. Transverse waves will tend to raise the cross 
correlation value at zero temporal offset wherever they are present and the emission from the 
strands which carry them is significant in both channels. This is simply a consequence of the 
wave bringing the strands (and their emission) in and out of the field of view of the pixel. 

Finally, we also observe a significant amount of zero time lag pixels in the 171-131 maps. This 
could indicate that though the plasma is cooling, it is not cooling much below the 171 peak 
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temperature of 0.8 MK. Alternatively, it may be that the plasma cools so rapidly from 0.8 to 0.5 
MK (for example, if enthalpy cooling dominates) that we cannot detect a time lag with the 30s 
resolution data that we use. A third possibility is that the contributions to the 131 channel near 
0.7 MK are much more than the current response functions indicate (e.g., Testa et al. 2011). 

Conclusion 

We measure time lags between the AIA coronal channels across the entire active region. 
These measurements provide crucial constraints on any model which seeks to describe the 
heating and subsequent cooling of coronal plasma. Our main findings are: 

1) We observe cooling plasma throughout the active region, for the entire 24 hour duration 
analyzed, independent of time series window size. The cooling plasma is observed both in 
discernable ‘loops’ as well as in the diffuse emission of the active region. Though we apply the 
same analysis to all of the pixels in the image without regard to AR feature, the basic AR shape 
is apparent in all of the maps. 

2) Time lags are observed between almost all channel pairs, including those with large 
temperature separation (e.g., 335-171), and time lags increase between channels with wider 
temperature separations. This indicates that coronal strands cool fully from temperatures greater 
than at least 3MK, and sometimes as high as 7 MK, down to temperatures lower than ~ 0.8 MK. 

3) The 94 results exhibit systematic behavior as a function of location in the active region. 
We find that the hot, 7 MK component of the 94 channel dominates the light curves of the core 
of the active region, while the cool, 1 MK component dominates outside the core, which suggests 
that nanoflare energies decrease away from the active region core (Viall & Klimchuk 2011). 
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4) Time lags are small for lines of sight that include moss, the transition region footpoints of 
hot (>3MK) coronal strands, which is also consistent with impulsive nanoflare heating. 

A significant amount of the emitting plasma is not steady; rather, it is dynamic and 
evolving. We discuss these results in the context of impulsive nanoflare heat release where the 
plasma cools fully before it is subsequently reheated, as shown in Viall & Klimchuk (2011). Our 
results are not consistent with the majority of the emission being caused by truly steady heating 
which produces steady emission. In a companion paper we examine these results in the context 
of effects of line of sight integration where physically separate strands along the same line of 
sight undergo out-of-phase impulsive heating. Finally, in this paper we have only analyzed a 
single active region, and there is no guarantee that other active regions, or even this active region 
later in its evolution, behave the same way. We are currently investigating other active regions in 
various phases of their evolution to address these questions. 
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Figures Captions 

Figure 1. Normalized temperature response functions (adapted from Viall & Klimchuk 2011) for 
6 SDO/AIAEUV channels: 131 (black), 171 (cyan), 193 (orange), 211 (blue), 335 (green) and 
94 (red). 

Figure 2. NOAA AR 1 1082 in six SDO/AIA channels at 3.5 UT on 19 June 2010. Upper left 
comer indicates channel and peak temperature sensitivity of that channel. From left to right, top 
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to bottom the images show channels 94, 335, 211, 193, 171 and 131. The images are displayed 
on a linear scale, and 1 pixel = 0.6”. Arrow indicates pixel used for Figure 3. 

Figure 3. Normalized light curves over two hours of location indicated with arrow in Figure 2 
images. From bottom to top we plot channels 131, 171, 193, 211, 335, and 94, each offset by 0.5 
in y. 

Figure 4. Cross correlation values as a function of offset time between pairs of light curves 
shown in Figure 3. 211-193 (green), 335-211 (blue); 335-193 (orange); 335-171 (cyan); 171- 
131(black) 94-335(red). Dots indicate time lags. 

Figure 5a. Peak cross correlation time lag maps for 0-12 UT, 19 June 2010 for the field of view 
shown in Figure 2. The color bar on the bottom indicates the time lag range in seconds. The 
channel pair is indicated on the top of each panel 

Figure 5b. Same as 5a. Note that the 211-193 and 171-131 pairs have different color bars. 

Figure 6a Peak cross correlation time lag maps for 0-2 UT, 19 June 2010. The color bar on the 
bottom indicates the time lag range in seconds. The channel pair is indicated on the top of each 
panel. 

Figure 6b. Same as 6a. Note the 211-193 and 171-131 pairs have different color bars. 
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Figure 7a. 335-211 time lag maps for all 2-hr time series. Temporal window in hours UT 
indicated on top. Top left panel is the same as 335-211 map in Figure 6a. 

Figure 7b Same as 7a. 
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